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Abstract 

In the framework of a multi-phase transport model (AMPT) with both partonic and hadronic 
interactions, azimuthal correlations between trigger particles and associated scattering particles 
have been studied by the mixing-event technique. The momentum ranges of these particles are 
3 < p%'^ < 6 GeV/c and 0.15 < ^5?^°^ < 3 GeV/c (soft), or 2.5 < p^*^ < 4 GeV/c and 
1 < ^5?^°"^ < 2.5 GeV/c (hard) in Au + Au collisions at ^/snn = 200 GeV. A Mach-hke structure 
has been observed in correlation functions for central collisions. By comparing scenarios with 
and without parton cascade and hadronic rescattering, we show that both partonic and hadronic 
dynamical mechanisms contribute to the Mach-like structure of the associated particle azimuthal 
correlations. The contribution of hadronic dynamical process can not be ignored in the emergence 
of Mach-like correlations of the soft scattered associated hadrons. However, hadronic rescattering 
alone cannot reproduce experimental amplitude of Mach-like cone on away-side, and the parton 
cascade process is essential to describe experimental amplitude of Mach-like cone on away-side. In 
addition, both the associated multiplicity and the sum of pT decrease, whileas the (px) increases, 
with the impact parameter in the AMPT model including partonic dynamics from string melting 
scenario. 



Key words: di-hadron azimuthal correlation, Mach cone, parton cascade, hadronic rescattering, 
AMPT 

PACS: 25.75.-q, 24.10.Nz, 24.10.Pa, 25.75.Ld 



I. Introduction 

A phase transition between hadronic mat- 
ter and the quark-gluon plasma (QGP) at 
a critical energy density of ~ 1 GeV/ fm^ 
has been predicted by Quantum Chro- 
modynamics (QCD) [1], which has moti- 
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vated the scientific program at the Rela- 
tivistic Heavy-Ion Collider (RHIC) at the 
Brookhaven National Laboratory. A very 
dense partonic matter has been shown to 
be produced in the early stage of central 
Au-hAu collisions at .Jsf^ = 200 GeV [2]. 
Many interesting phenomena have been 
observed including the measurements of 
elliptic flow [3], strangeness [4], J/-;/; [5] and 
jet quenching [6]. 
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Jet has been proved particularly good 
probe in RHIC experiments [6]. The high 
transverse momentum (px) partons (jets) 
which emerge from hard scattering pro- 
cesses will lose energy when they pass 
through the dense QCD medium. The en- 
ergy loss (jet quenching) mechanism results 
in distinct experimental observations such 
as the disappearance of one jet in back-to- 
back jet correlation at high [7]. At the 
same time, the loss energy must be redis- 
tributed in the soft pr region [8,9,10,11]. 
Experimentally the soft scattered parti- 
cles which carry the lost energy have been 
observed statistically via two-particle an- 
gular correlation of charged particles [12]. 
Reconstruction of these particles will con- 
strain models which describe production 
mechanisms of high px particles, and may 
shed light on the underlying energy loss 
mechanisms and the degree of equilibra- 
tion of jet products from energy loss in the 
medium. 



A Mach-likc structure (the splitting of the 
away side peak in di-jct correlation) has 
recently been observed in azimuthal corre- 
lations of scattered secondaries associated 
with the high Pt hadrons in central Au+Au 
collisions at ^/s^ = 200 GeV [13,14,15]. 
Several theoretical interpretations have 
been proposed for this phenomenon. For 
instances, Stocker proposed Mach cone 
structure from jets traversing the dense 
medium as a probe of the equation of state 
and the speed of sound of the medium [16]. 
Casalderrey-Solana, Shuryak and Teaney 
argued for a shock wave generation because 
jets travel faster than the sound in the 
medium [17]. They fit the broad structure 
on the away side of the azimuthal correla- 
tion with a Mach cone (shock wave) mecha- 
nism; Koch and Wang et al. could produce 
a Mach-like structure with a Cherenkov 
radiation model [18]; Armesto et al. inter- 
preted the sideward peaks result of 
medium dragging effect in Ref. [19]; Rup- 
pert and Miiller argued that the Mach-like 
structure can appear due to the excitation 
of collective plasmon waves by the moving 
color charge associated with the leading 
jet [20]; Renk and Ruppert applied a re- 
alistic model for the medium evolution to 



explain the observed splitting of the away 
side peak in Ref. [21]; Chaudhuri studied 
the effect on Mach-like structure from jet 
quenching on the hydrodynamical evolu- 
tion of a QGP fluid [22]; Satarov, Stocker 
and Mishustin investigated Mach shocks 
induced by partonic jets in expanding 
quark-gluon plasma [23]; A conical flow 
induced by heavy-quark jets was also pro- 
posed by Antinori and Shuryak [24]; Hwa 
et al. discussed the centrality dependence 
of associated particle distribution for Au 
+ Au and d + Au, respectively, in recom- 
bination model [25] . The dynamical nature 
of the Mach-hke structure continues to be 
a subject of many theoretical and experi- 
mental investigations. 

There is no quantitative interpretation 
from dynamical transport models yet for 
the Mach-like structure in particle corre- 
lations from hard scattering particles in- 
teracting with the dense medium. In this 
work, we report a study of associated par- 
ticle correlations with a triggered particle 
and investigate the Mach-like structure 
using A Multi-Phase Transport model( 
AMPT ) [26] . We have applied the mixing- 
event technique in the AMPT analysis 
as what has been used in the analysis of 
RHIC collision data to remove background. 
We found that both parton cascade and 
hadronic rescattering can produce appar- 
ent correlations between triggered and as- 
sociated particles similar to the Mach-like 
structure. But the hadronic rescattering 
mechanism alone is not able to produce 
a large enough amplitude for Mach-like 
cone on away side, and the parton cascade 
process is indispensable. 

The paper is organized as follows. In Sec- 
tion II, we give a brief description of the 
AMPT model and the initial conditions. 
Section III describes mixing-event tech- 
nique in our simulation analysis. Results 
and discussions are presented in Section 
IV. Finally a summary is given in Section 
V. 

II. Brief description of AMPT Model 



The AMPT model [26] is a hybrid model 
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which consists of four main components: 
the initial conditions, partonic interac- 
tions, the conversion from partonic matter 
into hadronic matter and hadronic rescat- 
tering interactions. The initial conditions, 
which include the spatial and momentum 
distributions of minijet partons and soft 
string excitation, are obtained from the HI- 
JING model [27]. Excitation of strings will 
melt the string into partons. Scatterings 
among partons are modelled by Zhang's 
parton cascade model (ZPC) [28], which 
includes two-body scatterings with cross 
sections obtained from the pQCD calcula- 
tions with screening mass. In the default 
AMPT model [29] partons are recombined 
with their parent strings when they stop 
interacting, and the resulting strings are 
converted to hadrons by using the Lund 
string fragmentation model [30]. In the 
AMPT model with the option of string 
melting [31], a quark coalescence model is 
used to combine partons into hadrons. Dy- 
namics of the subsequent hadronic matter 
is then described by A Relativistic Trans- 
port (ART) model [32]. Details of the 
AMPT model can be found in a recent re- 
view [26]. Previous studies [31] have shown 
that the partonic effect could not be ne- 
glected and the string melting AMPT is 
much more appropriate than the default 
AMPT version when the energy density is 
much higher than the critical density for 
the predicted phase transition [26,31,33]. 
In the present work, the parton interaction 
cross section in AMPT model with string 
melting is assumed to be lOmb. 

III. Analysis Method 

In order to compare with experimental 
measurements of correlations between 
trigger and associated hadrons, we use 
the mixing-event technique to subtract 
combinatorial background in our analy- 
sis. Two ranges of pt window selections 
for the trigger and associated particles 
have been used: one is 3 < p^'^ < 6 
GeV/c and 0.15 < < 3 GeV/c 

(we call this selection as "soft" associ- 
ated hadrons since soft particles dominate 
in this Pt range for associated particles); 

the other is 2.5 < p^*^ < 4 GeV/c and 
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Fig. 1. (a): The A(j) distribution between a 
triggered hadron (3 < p^^^^ < 6 GeV/c) and 
associated hadrons (0.15 < p^'"'°'' < 3 GeV/c) 
(circles), where the background (dash line) has 
not been subtracted. The data are 200 GeV 
Au + Au collisions with 20-40% collision cen- 
trality from the AMPT model. The dash area 
is the region of ZYAM normalization (see texts 
for detail); (b): The A(j) distribution where 
the background has been subtracted by a mix- 
ing-event technique, where the solid line is a 
two-Gaussian fit. 

1.0 < p^''"" < 2.5 GeV/c (denoted as 
"hard" associated hadrons since there are 
more hard particle components than that 
in the "soft" selection). The triggered par- 
ticles and associated particles ("soft" and 
"hard" ones) both are selected with a 
pseudo-rapidity window |?7| < 1.0. In the 
same events, pairs of associated particles 
with a triggered particle are accumulated 
to obtain A0 = (f> — (ptrig distributions. In 
order to construct the background which 
is mainly from the effect of elliptic flow 
[12,14], a mixing-event method is ap- 
plied to simulate the background. In this 
method, we mixed two events which have 
very close centrality into a new event, and 
extracted A0 distribution to be used as 
background distribution. When subtract- 
ing the background from the same events, 
ZYAM (zero yield at minimum) assump- 
tion is adopted which has been used in the 
experimental analysis [14]. Figure 1 shows 
the A(j) distribution for trigger hadrons 

of 3 < p^*^ < 6 GeV/c and associated 
hadrons of 0.15 < < 3 GeV/c before 

and after the background subtraction. The 
data are from AMPT model simulation of 
Au -|- Au collisions at ^/sNN = 200 GeV 
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Fig. 2. Soft scattered associated hadron 
A(f) correlations for trigger hadrons of 
3.0 < p^*^ < 6.0GeV/c and associated 
hadrons of 0.15 < < 3.0 GeV/c 

from AMPT Monte Carlo simulations of 
Au + Au collisions at ^/snn = 200 GeV 
with various collision centralities. Triangles: 
string melting version after hadronic rescat- 
tering; circles: string melting version before 
hadronic rescattering; diamonds: default ver- 
sion after hadronic rescattering; stars: de- 
fault version before hadronic rescattering; 
squares: experimental data are from Ref [12] 

where 4.0 < p^*^ < 6.0 GeV/c and 
0.15 < < 4.0 GeV/c. 



In Ref. [12], it was found that the A0 = 
— (ptrig distribution of recoiling hadrons 
from a high px triggered particle is signifi- 
cantly broadened in central Au -|- Au colli- 
sions at y/SNN = 200 GeV, which supports 
the picture of the dissipation of jet energy 
in the medium. In order to increase the sta- 
tistical sample of triggered particles in our 
calculation, we set px range for trigger par- 
ticles to 3 < p^*^ < 6 GeV/c and for asso- 
ciated particles to 0.15 < p^''°'' < 3 GeV/c 
in our analysis. Both triggered and associ- 
ated particles are further selected with a 
pseudo-rapidity cut of \ri\ < 1.0. 

Figure 2 presents AMPT model calcula- 
tions of soft associated hadron Acp corre- 
lations from Au + Au collisions at ^/sNN 
= 200 GeV for different centralities un- 
der various AMPT running conditions. 
In order to compare our results with ex- 
perimental data which give the correla- 
tions among associated charged hadrons, 
the experimental data are multiplied by 
a factor of 1.58 to account for the contri- 
bution from neutral hadrons [12,35]. Note 
that for the default AMPT version before 
hadronic rescattering, we here give A(j) cor- 
relation for 0-10% collision centrality only 
due to the lack of statistics obtained for 
other centrality bins. We found that the 
hadronic rescattering increases A0 correla- 
tion yields for both versions. A very strong 
Mach-like structure is observed for central 
Au-|-Au collisions before hadronic rescat- 
tering in the melting AMPT version, which 
indicates that the Mach-like structure has 
been formed in parton cascade process. 
The effect of hadronic rescattering on the 
A(j) correlation in central collisions (< 20% 
centrality) does not wash out, even slightly 
enhances, the Mach-like structure, which 
is qualitatively in agreement with the ef- 
fect of time-dependent speed of sound on 
the development of the conical wave in 
expanding QCD matter [36]. 



for 20-40% collision centrality, where the 
string melting mechanism has been turned 
on. 

IV. Results and Discussions 



Figure 3 shows the AMPT calculations of 
A0 correlations between a triggered par- 
ticle of 2.5 < p^"^ < 4.0GeV/c and hard 
associated particles of 1.0 < p^^"""" < 2.5 
GeV/c from Au + Au collisions at ^/snn 
= 200 GeV for various centralities, where 
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Fig. 3. AMPT model calculations of A(/> 
correlations between a triggered hadron of 
2.5 < p^*^ < 4.0 GeV/c and hard associated 
hadrons of 1.0 < ■p9^^°'^ < 2.5 GeV/c from Au 
+ Au collisions at y/sJ/N = 200 GeV in various 
centralities. Triangles: AMPT version with 
string melting after hadronic rescattering; Full 
circles: AMPT version with string melting be- 
fore hadronic rescattering; Diamonds: default 
AMPT version after hadronic rescattering; 
Stars: default AMPT version before hadronic 
rescattering; Open circles: experimental data 
from Ref . [14] ; Hatched areas indicate the ex- 
perimental uncertainty. Note that the scaling 
factors may also reflect different pseudorapid- 
ity ranges used in experimental data analysis 
and model calculations. See text for details. 

Dseudorapidity cuts of < 1.0 and 

^assoc^ < 1.0) have also been applied. We 
bund that the effect on A0 correlations 
from hadronic rescattering is much smaller 
than the case for soft associated particles, 



which may indicate that fewer fraction of 
hard associated hadrons suffer hadronic 
rescattering. In addition, Mach-like struc- 
tures have been observed on away-side cor- 
relations in both string-melting and default 
AMPT versions. In the default AMPT ver- 
sion, however, the Mach-like structures can 
be observed only after turning on hadronic 
rescatterings. 

The default AMPT version appears to pro- 
duce the number of associated particles 
matching better with the experimental 
measurement than that from the string- 
melting AMPT version. This may be at- 
tributed to the fact that the string-melting 
AMPT model always produces softer pt 
spectra than the default AMPT model be- 
cause current quark masses have been used 
in the partonic cascade stage [26]. It may 
be improved if thermal parton masses are 
applied [34] . However, in the present study, 
we only focus on the correlation shape in- 
stead of yields in the correlation region. In 
this context, we shall compare the shapes 
of A0 correlation functions between the 
AMPT model and experimental data. To 
this end, the AMPT model data have been 
multiplied with different normalization fac- 
tors which are listed on right sides of each 
panel in Figure 3. Please note that there 
is a difference for pseudo-rapidity cut be- 
tween our simulation and the experimental 
data. Our simulation cut of associated par- 
ticle is \r]\ < 1.0 (due to the limited statis- 
tics), while \t]\ < 0.35 in PHENIX data. In 
this case, the scaling factors of our simula- 
tion should be reduced by a factor of about 
3 in order to match the data. By this reduc- 
tion, the scaling parameters shown in the 
figures have been normalized to the same 
\ri\ cuts as the data. However, even with 
the scaling of 3 these factors still deviate 
from 1 and show a centrality dependence. 
These deviations can be attributed to the 
following facts: (1) a large yield of particles 
in the correlation region can partially come 
from the fact that our pt spectra in AMPT 
calculation with string melting scenario 
is under-predicted in high p^ region if we 
compare with the experimental data. In 
this case, the soft component will be over- 
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Fig. 4. Splitting parameter D versus im- 
pact parameter for Au + Au collisions 
at y/SNN = 200 GeV {pT windows in 
model: 2.5 < pp^ < 4.0 GeV/c and 
1.0 < p^"""" < 3 GeV/c; pr windows in 
experiments: 2.5 < p^*^ < 4 GeV/c and 
2.0 < < 3.0 GeV/c). Open symbols: the 

string-melting version with hadron rescatter- 
ings; half filled symbols: the default version 
with hadronic rescatterings; filled symbols: ex- 
perimental data from Ref. [15]. 

the data and simulation; (2) the centrality 
dependence of scaling factors may stem 
from more excessive parton interactions 
when a jet parton passes through various 
dense partonic matter depending on the 
collision centrality. In this study, our em- 
phasis will be on the comparison of shapes 
of the Mach-like structure between the 
AMTP model calculations and experimen- 
tal data. Our simulation results indicate 
that the string-melting AMPT model can 
describe shapes of A0 correlations between 
a triggered hadron and associated parti- 
cles better than the default AMPT model, 
especially for central Au+Au collisions. 

In order to quantitatively characterize the 
Mach-like structure, a splitting parameter 
D has been extracted in our analysis. The 
splitting parameter D is defined as half 
distance between two Gaussian peaks on 
away side of associated particle A0 cor- 
relations. The D value reflects the size of 
direction-splitting of Mach-cone on the 
away side. Figure 4 shows the impact pa- 
rameter dependence of D. Our results in- 



dicate that the string-melting AMPT ver- 
sion can roughly match the experimental 
data. Note that of our px cuts are slightly 
different from experimental cuts in order 
to increase our statistics. However, D val- 
ues from the default AMPT version are 
significantly smaller than the experimental 
data. We conclude that hadronic rescat- 
tering mechanism alone is not enough to 
produce the amplitude of Mach-like cone 
structure on the away side and parton 
cascade mechanism is necessary. It should 
be noted that minijet partons produced 
from hard scatterings can lose energy by 
gluon radiations and transfer their ener- 
gies to nearby soft strings in the HIJING 
model, which has been called jet quenching 
[27]. In the AMPT model, the jet quench- 
ing in the HIJING model is replaced by 
parton scatterings from the ZPC. Only 
two-body scatterings are included in the 
current ZPC, and higher-order contribu- 
tions such as 2^3 scattering processes 
to the jet energy loss are still missing in 
our AMPT calculations. It has been pro- 
posed that higher order processes may 
contribute significantly to the observed 
large jet quenching [37]. However our re- 
sults indicate that 2^2 processes could 
account for most of the amplitude for the 
Mach-like structure. The AMPT version 
with string melting scenario has also been 
shown to generate elliptic flow of hadrons 
better than the default version. Further- 
more, it can reproduce the mass ordering 
in the particle dependence of elliptic flow, 
which can be described by hydrodynamics 
models [31,33]. It is believed that the big 
cross section for parton interactions used 
in the AMPT model leads to strong par- 
ton cascades that couple partons together, 
resulting the onset of hydrodynamics be- 
havior from parton cascades [38]. It has 
also been proposed that a Mach-like struc- 
ture can be generated if the jet velocity is 
faster than the sound velocity when the jet 
traverses through the dense matter [16,17]. 
In addition, continuous partonic rescatter- 
ings and the rapid expansion of the dense 
medium can also lead to parton's direction 
in the medium to both sides of triggered 
particle direction, which may result in an 
apparent Mach-like structure in the AMPT 
model. Nevertheless, we conclude that sig- 
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nificant parton cascades are essential in 
order to describe elliptic flow of liadrons 
and the Mach-like structure in particle A0 
correlations simultaneously. 

We have also investigated several charac- 
teristics of soft associated particles in the 
AMPT model with the string-melting sce- 
nario. Figure 5 shows a comparison be- 
tween experimental data and our AMPT 
calculation results with hadronic rescatter- 
ings. Panel (a) in Fig. 5 shows the num- 
ber of associated hadrons function of 
impact parameter in Au+Au collisions at 
^/sNN = 200 GeV. The trigger particles 
are not included in the Nch counting for 
the near side. The number of average as- 
sociated hadrons decreases with increasing 
impact parameters in the model as well as 
in the experiment. However, the number of 
associated hadrons in the AMPT model is 
much larger than the experimental data on 
both near and away sides, which may due 
to different pt cuts used in the model cal- 
culation and experimental data. Panel (b) 
shows the pt distributions of the associated 
hadrons on near side and away side in most 
central (0-5% centrality) Au + Au colli- 
sions at ^/sNN = 200 GeV. Panel (c) gives 
the sum of pt magnitude, which approxi- 
mates to associated energy, as a function 
of impact parameter. The triggered parti- 
cles are included in the pt magnitude sum 
for the near side. The sums of pt magni- 
tude decrease with impact parameters for 
both near side and away side, and theoreti- 
cal values are somewhat higher than the ex- 
perimental data, again may be affected by 
the different pt cuts. Panel (d) displays the 
dependence of (pt) on impact parameters 
on away side, (pt) increases with impact 
parameter. These comparisons may help to 
address the issue of parton thermalization 
through the parton cascade mechanism in 
central Au+Au collisions [12] . 

V. Conclusions 

In summary, the origin of the Mach-like 
structure in correlations between triggered 
hadrons and soft or hard associated par- 
ticles has been investigated in the frame- 
work of a hybrid dynamics transport model 
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Fig. 5. Comparisons between model calcula- 
tions {pt windows: 3 < p^*^ < 6 GeV/c and 
0.15 < p^^°'^ < 3 GeV/c) and experimental 
data [12] {pT windows: 4 < p^*^ < 6 GeV/c 
and 0.15 < pf'"'' < 4 GeV/c) in Au + Au col- 
lisions at y^sjvTV = 200 GeV. Open symbols: 
calculations from the string-melting AMPT 
version, filled symbols: experimental data, cir- 
cle: near side, triangle: away side. Panel (a): 
Number of associated hadrons versus impact 
parameter; Panel (b): pT distributions for the 
associated hadrons on near side and away side 
for 0-5% centrality; Panel (c): Sum of mag- 
nitude as a function of impact parameter on 
near side and away side; Panel (d): (px) as a 
function of impact parameter on away side. 

which includes two dynamical processes, 
namely parton cascade and hadronic 
rescattering. By comparing the different 
calculation results with or without parton 
cascade, before or after hadronic rescatter- 
ing, we found that the associated particle 
correlations and the Mach-like structure 
have been formed mostly before hadronic 
rescatterings, which indicates that these 
kinds of correlations are born in the par- 
tonic process and are further developed in 
the later hadronic rescattering processes. 
For the hard associated particles, hadronic 
rescattering hardly changes the Mach-like 
structure which is mostly formed in par- 
ton cascade processes. Therefore, charac- 
teristics of hard associated particles may 
directly reflect information of intrinsic par- 
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tonic dynamics. Meanwhile, the effect of 
hadronic rescattering can not be ignored 
especially for soft associated particles, be- 
cause many of these soft associated parti- 
cles either are produced in or suffer from 
hadronic rescattering processes. Our calcu- 
lations indicate that hadronic rescattering 
mechanism alone is unable to produce a 
splitting parameters D distribution for the 
Mach-like structure matching the experi- 
mental data. The parton cascade mecha- 
nism is essential for the Mach-like structure 
while the exact shape and strength of the 
azimuthal correlations seem to depend on 
detailed properties of the partonic medium 
and jet, which still awaits quantitative 
explanations. 
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